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Abstract 

The tensile yielding of multiwall carbon nanotubes (MWCNTs) has been 
studied using Molecular Dynamics simulations and a Transition State Theory 
based model. We find a strong dependence of the yielding on the strain rate. 

A critical strain rate has been predicted above/below which yielding strain of 
a MWCNT is larger/smaller than that of the corresponding single- wall carbon 
nanotubes. At experimentally feasible strain rate of 1% /hour and T = 300K, 
the yield strain of a MWCNT is estimated to be about 3-4 % higher than 
that of an equivalent SWCNT, in good agreement with recent experimental 
observations. 
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Carbon Nanotubes (CNTs) have been subject of extensive studies in recent years. In 
experiments CNTs are grown in arc-discharge, chemical vapor deposition, and laser ablation 
based methods. [1] Singlewall CNTs (SWCNTs) and Multiwall CNTs (MWCNTs) with 
many concentric shells of graphene type C layers are the common products depending on 
experimental conditions. Theoretical and experimental studies have shown that CNTs have 
exceptionally strong mechanical characteristics, wuth Young’s modulus as high as 1 TPa. 
Many applications have been suggested to exploit the mechanical properties of CNTs and 
MWCNTs have been the main candidate as they can be produced in large quantities in 
experiments. The understanding of the mechanical response, such as yielding behavior, of 
MWCNTs is thus significantly important from the applications perspective. 

Based on extensive Molecular Dynamics (MD) simulations of yielding of SWCNTs, re- 
cently we have proposed a transition state theory (TST) based model for the tensile failure 
of SWCNTs. [2] In the model the tensile yielding strain is considered to depend on both 
the temperature and the strain rate. Starting with an Arrhenius description, the activation 
time for a yielding process is written as t = ±e E "/ kBT , where E u is the activation energy 
and v is vibrational frequency. This can be expanded to t = k e {Ev-VKt)/k B T because E„ 
linearly decreases with applied tensile strain e with K as the force constant and V as the 
activation volume. For the case of SWCNT the tensile yielding is triggered by Stone- Wales 
bond rotation process, and E v is the activation energy of formation such defects (References 
on SW rotation and its static activation energy of formation are given in [3].) Detailed 
analysis shows that the yielding strain ey can be obtained by inverting the above expression 
for the activation time as according to, [2] 
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where E„ is the averaged dynamic activation energy of SW bond rotations; eo is the 
intrinsic strain rate; n site is the number of atomic sites available for defect formation; and 
N is the number of SW bond rotations involved in the yielding process. The yield strain 
of a 60 A long CNT(10,0) as a function of strain rate at different temperatures, shown in 
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Fig. 1, was computed from extensive MD simulations (to be discussed below), and found 
to follow the relation expressed in Eq. (1). The intrinsic parameters obtained from the MD 
simulations data [2] for a CNT(10,0) are as following: E v = 3.6eV; e 0 /N = 8 x 10~ 3 ps -1 ; 
and V = 2.88A 3 which is similar to that of the vacancy volume on a C atom on CNT, and 
will be shown to play an important role in the yielding of MW CNT as well. 

The dynamics and yielding of MWCNTs is the subject of interest in this paper. As a 
model, we consider a MWCNT made of two shells, CNT(20,0) (63.9 A, 1200 atoms) and 
(10,0) (55.2 A, 524 atoms) (Fig. 2). The capped inner shell is put concentrically and 
symmetrically inside the outer shell with inter shell spacing of 3.9 A. The tensile strain is 
applied at the two ends of the outer shell gradually (0.25% tensile strain per step followed by 
dynamics relaxation for a chosen time interval which is dependent on the strain rate) . Any 
inter shell load transfer is through VDW interactions [4]. Langevin friction force scheme 
[6] is used to control temperature and time step is kept to 0.5 femto second. The Tersoff- 
Brenner potential is used for C-C interactions, which can describe C bond formation and 
breaking. Details of this potential are described elsewhere [5]. 

The MD simulations are conducted at a high temperature of 2400K. Due to the tensile 
strain, the MWCNT (20,0)/(10,0) undergoes a yielding process initiated by the formation of 
SW T bond rotation defects followed by the breaking of the outer shell, while the structure of 
the inner shell remains intact. The atomic structures are monitored and recorded during the 
dynamics. The yield strain is defined as the tensile strain at which the MWCNT behaves 
differently from its elastic behavior with an abrupt drop of the strain energy. The yield 
strain of MWCNT (20,0)/(10,0) as a function of strain rate, in comparison with that of a 
SWCNT (20,0) under similar conditions, is shown in Fig 3. The standard error at a chosen 
point as shown in Fig. 3 is calculated from a data of 5 sample sets. It is seen that the 
MWCNT has the similar logarithmic dependence on the strain rate as was observed for 
SWCNTs in our earlier work. [2] This is not surprising because we have assumed that the 
entire tensile strain is applied on the single outermost shell. The role of the inner shells in 
a MWCNT is to modify this behavior through VDW forces. 
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One important feature from Fig. 3 is that there is a critical crossover point e cross of the 
strain rate c at which yield curves in the two cases cross. For c faster than Across? the MWCNT 
has a larger yield strain compared with the equivalent SWCNT. The situation is reversed 
for e slower than e cro3S . This feature can be explained by examining the intrinsic parameters 
that appear in Eq. (1). Experiments have suggested that for both dispersed SWCNT ropes 
and MWCNTs, the Young’s modulus is similar at 1 TPa. [1] The intrinsic strain rate and 
the number of processes involved in yielding are also expected to be of similar magnitude. 
A most probable parameter inducing the observed difference thus could be the activation 
volume because the magnitude of which affects the slope and the intersection of the yield 
strain vs strain rate plot in opposite ways. This results in a critical crossover point of the 
strain rate. By fitting the two curves in Fig. 3, we have (1) for a (20,0)/ (10,0) MWCNT, 

= 0.2339 + 0.0148 Inc; and (2) for a (20,0) SWCNT, ey = 0.1802 + 0.0094 In e. Assuming 
other parameters remaining same the ratio of the activation volume on a MWCNT to that 
on a SWCNT is found to be jjjgj = 0.64 or = 0.77 by comparing the slopes or the 
intersections, respectively. The two estimated values roughly agree with each other. The 
activation volume on the outermost shell of a MWCNT is thus found to be within 65 to 75% 
of that on an equivalent SWCNT shell. 

A physical reason for the reduction of the activation volume on the outermost shell of 
a MWCNT can be explained as following. In the formation of SW bond rotation, a C-C 
bond rotates around its center by 90 degree to form a heptagon-pentagon pair. It has been 
shown that the out of plane rotation of the bond during the process provides a lower barrier 
pathway as compared to that of a totally in-plane rotation [7]. The out of plane rotation 
of the C-C bond on the outer shell in a MWCNT is thus restricted to a smaller activation 
volume due to the VDW forces from the inner shells. 

The critical crossover point in general cases can be estimated as following. To have 
same yield strain for systems with different activation volume, assuming other parameters 
remaining same, the following equation must hold for the strain rate at the crossover point 

Across j 
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Where e CT0SS decays exponentially with temperature and is m the range of day 1 or 
smaller only for T above 1000K. At T = 300K, as shown in Fig. 3, e cross is in the range of 
10 -^ year -'' Thus at room temperature, for all practical purposes, the yield strain of the 
outermost shell of a MWCNT is always considered to be larger than that of a SWCNT. 

Above MD simulations were conducted at T = 2400K. In usual experiments, the tem- 
perature is at 300K, and strain rate is much slower, which is in the range of minute 1 or 
hour- 1 . There is an equivalence between high temperatures with fast strain rates and low 
temperatures with slow strain rates. By rewriting Eq. (1), such equivalence can be expressed 

as following, 

jiW_) Tl = ( (4) 
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Where Eq. (4) can be used to accelerate the yielding processes in the MD simulation 
by performing the simulations at high temperature and fast strain rates. In Fig. 3, the 
equivalent strain rate at T = 300K, as according to Eq. (4) and with parameters listed in 
second paragraph with n slte = 1200 for the MWCNT used in our simulation, is marked on 
the top axis of the figure. The range marked as an ellipse in Fig. 3 is for the strain rate 
of hour- 1 at T = 300K, which is usually used in experiments. It can be seen that in this 
range, the yield strain of the outermost shell of the MWCNT is about 3-4 % higher than that 
of the equivalent SWCNT shell. This is in qualitative agreement with recent experimental 
observations, where SWCNTs have been found to yield at up to 6% tensile strain [8,9] and 

MWCNTs yield at 10-12% tensile strain [10]. 

The change of the strain energy of the inner shell (10,0) is plotted as a function of the 
tensile strain on the outer shell (20,0) in Fig. 4 for four strain rate cases at T = 2400K. 
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For small strains within 2%, all the cases show the sign of load transfers. For strains larger 
than that, the load transfer begins to cease, which suggests that the VDW forces between 
the shells is not strong enough to hold additional interfacial tensile stress. The strain energy 
curves fluctuate (possibly due to thermal fluctuations in the simulation) more for some strain 
rates than the others in this range and no clear trend is observed. With the tensile strain 
continuing to increase to a yield point, there is a sudden increase of load transfer (shown as 
the abrupt increase of strain energy at higher strains in Fig. 4), indicating the necking of 
the out shell during yielding process, which causes an pressure onto the the inner shell and 
helps more load transfers. 

In summary, a recently proposed Transition State Theory based model of tensile Welding 
of SWCNTs has been extended to explain the yielding of MWCNTs with the assumption that 
due to contact all the tensile strain is applied only to the outermost shell of the MWCNT. The 
yield strain of a two shell (20,0) / (10,0) MWCNT are computed through MD simulations and 
are found to decrease logarithmicly with slower strain rate, in agreement with the underlying 
TST based model. A critical crossover strain rate is found below/above which yield strain 
of a MWCNT is smaller/larger than the equivalent SWCNT shell. This is attributed to 
the difference in the activation volume in the two cases, and the activation volume on the 
outermost shell of a MWCNT is found to be about 65 to 75 % of that on an equivalent 
SWCNT. At room temperature, and experimentally feasible strain rates, the tensile yield 
strain of a MWCNT is found to be 3-4 % larger than the equivalent SWCNT which is in 
good agreement with recent experimental observations. 
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FIGURES 

FIG. 1. The Yield strain of a 60 A long CNT(10,0) as a function of strain rate at various 
temperatures from 300K to 2400K from MD simulations. 

FIG. 2. The two-shell MWCNT consisting of CNT (20,0) and (10,0). The inner shell is capped 
at its two end. 

FIG. 3. The tensile yield strain of the (20,0)/ ( 10,0) MWCNT (solid circle and line) as a function 
of strain rate at T = 2400K (bottom axis), in comparison to that of a CNT(20,0) (diamond and 
dashed line). The error bar at the point of strain rate of 2.5 x 10 - V? -1 i s calculated from data 
of 5 sample sets. The arrow in the lower left corner shows the crossover point. The strain rate in 
the top axis is corresponding to T = 300K, on which the marked positions shows rates of 1 %/hour 
and 1 Vo /year. The dotted ellipse shows the range of experimental feasible conditions with strain 
rate as 1 %/hour at T = 300K. 

FIG. 4. The change of the strain energy of the inner shell CNT(10,0) as a function of the 
applied tensile strain on the outer shell CNT(20,0) at four strain rates at T=2400K. 


8 



Strain rate (1/ps) 









Tensile strain on outer shell (20,0) (%) 



2400K 




